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ABSTRACT 


The  NATO  Reference  Mobility  Model  (NRMM)  produces  a  prediction 
of  the  maximum  average  speed  at  which  a  vehicle  can  traverse  an  area 
(terrain  unit) .  The  program  described  here  uses  selected  values 
calculated  in  the  NRMM  to  determine  the  factor  which  is  the  limiter 
of  speed  for  a  vehicle  and  terrain  unit.  In  the  case  of  a  NO-GO 
prediction,  the  reason  for  the  NO-GO  is  deduced.  Detailed  and  summary 
diagnostic  tables  are  produced  together  with  a  graphical  presentation 
of  the  diagnostics. 
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INTRODUCTION  AND  BACKGROUND 


The  NATO  Reference  Mobility  Model  (NRMM)  [1]  is  a  computer 
program  which  provides  a  comprehensive  assessment  of  cross-country 
performance  of  a  vehicle.  The  basic  output  of  the  NRMM  is  the  speed 
at  which  a  vehicle  can  travel  in  a  terrain  region  (patch)  judged 
uniform  with  respect  to  mobility.  (This  speed  is  often  referred 
to  as  "Speed-Made-Good"  since  the  path  need  not  be  a  straight  line 
due  to  mar  auvers  around  trees  and  obstacles . )  This  speed-made-good 
for  a  patch  is  calculated  as  an  interaction  of  a  variety  of  descriptors 
of  the  terrain,  vehicle  and  scenario  and  includes  assessment  of  the 
effects  of  such  things  as  the  vehicle  powertrain  and  geometry,  the 
strength  of  the  soil  in  the  patch,  vegetation,  driver  dictated  limits, 
esc.  In  addition  to  the  single  output  number  of  speed-made-good, 
results  of  literally  hundreds  of  intermediate  computations  may  be 
obtained  from  the  NRMM. 

The  objective  of  this  study  was  to  identify  those  intermediate 
computations  in  the  NRMM  which  could  be  used  to  give  to  the  user 
better  insight  into  the  vehicle/terrain  interaction  without  the  need 
for  total  immersion  in  a  sea  of  numbers.  The  identification  of  an 
intermediate  result  as  useful  for  insight  has  been  governed  by  prior 
use  of  the  NRMM  and  the  needs  expressed  by  the  users  during  and  after 
these  prior  studies. 

In  particular,  the  way  in  which  the  model  has  been  used  (for 
vehicle  studies)  in  the  last  few  years  is  the  following: 

1.  Several  regions  in  the  world  have  been  identified  as 

being  of  interest  to  the  military  user  community  for 
mobility  study  purposes. 

2.  Available  terrain  data  has  been  analyzed  to  yield  a 

description  of  these  regions  as  a  mosaic  of  patches 
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(areas  judged  uniform  with  respect  to  mobility)  overlaid 
with  roads,  trails,  rivers  and  other  features  which 
appear  as  curves  and  lines  on  a  map. 

3.  The  values  of  the  terrain  descriptors  (soil  strength, 

topographic  slope,  roughness,  etc.)  have  been  established 
and  amassed  into  computer  data  files. 

4.  The  NRMM  was  run  to  obtain  a  speed-made-good  for  each 

terrain  unit  (patch,  road,  or  trail)  in  these  files. 

5.  This  set  of  speeds  was  processed  to  obtain  a  few  performance 

measures  —  usually  percentage  of  terrain  denied  (zero 
speed  or  NO-GO)  and  the  average  speed  (weighted  by  area) 
on  the  least  severe  90%  of  the  terrain  (without  regard 
to  the  spatial  distribution  of  the  most  severe  terrain) . 

A  speed  profile  is  often  presented  as  well. 

While  these  data  are  useful  in  comparing  vehicles,  in  certain 
studies  (WHEELS  [2],  HIMO  [3])  which  used  predecessors  of  the  NRMM, 
information  on  the  causes  of  NO-GO's  and/or  the  factors  which  limit 
the  speed  were  judged  to  be  required  and  were  produced.  The  program 
developed  in  this  study  produces  similar  diagnostic  output  for  runs 
of  the  NRMM  and  was  developed  with  some  attention  being  paid  to 
program  portability. 

In  the  NRMM,  various  single-factor  aspects  of  the  terrain- 
vehicle  interaction  which  can  cause  a  NO-GO  condition  or  limit  speed 
are  evaluated  (e.g.,  limitation  due  to  driver  tolerance  or  lack  of 
braking  ability)  and  then  multi-factor  combinations  of  effects  are 
assessed.  These  calculations  are  performed  (when  appropriate) 
traveling  with  and  against  the  input  topographic  slope  (up  and  down) 
and  (in  the  Areal  Module)  at  zero  topographic  slope  (level).  When 
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one  of  the  calculations  produces  a  NO-GO,  many  subsequent  calculations 
are  skipped;  otherwise,  a  variety  of  candidate  speeds  which  reflect 
the  effects  of  some  of  the  constraints  are  computed  and  then  the 
feasibility  speed  is  obtained.  If  any  of  the  up,  level  and  down 
speeds  is  zero,  the  speed-made-good  for  the  terrain  unit  is  set  to 
zero,  otherwise  the  harmonic  average  is  the  output. 

RESULTS  OF  CURRENT  STUDY 


Since  the  factor  which  controls  the  speed  may  be  different 
when  traveling  with  the  slope,  on  the  level  and  against  the  slope, 
it  was  decided,  for  this  study,  to  analyze  these  cases  separately. 

Another  consideration  in  setting  up  this  analysis  was  that  the 
mechanism  for  making  changes  to  the  NRMM  is  time-consuming,  so  that 
all  modifications  required  for  the  output  analysis  were  restricted 
to  the  Control  and  Input/Output  (C&I/O)  module  in  which  local  changes 
are  permitted  by  the  rules  governing  the  use  of  NRMM  established 
by  the  NATO  Technical  Management  Committee.  The  organization  of 
the  Areal  Module  of  the  NRMM  (with  its  own  control  subprogram)  allows 
the  changes  required  to  be  minimal  and  concentrated  in  a  few  subroutines 
of  the  C&I/O  Module.  These  changes  are  fully  described  later.  Creation 
of  the  same  type  of  output  analysis  for  the  Road  Module  is  straightforward 
but  requires  changes  in  the  Module  itself  since  there  is  no  separate 
control  subprogram  for  this  Module. 

The  output  analysis  is  performed  as  follows: 

1.  The  additions  to  the  NRMM  Control  and  Input/Output  Module 

are  entered.  (This,  of  course,  is  done  only  once.) 

2.  The  NRMM  is  run  for  the  desired  vehicle  and  terrain  with 

the  new  control  variable  KDIAG  set  to  1  which  produces 
the  additional  output  required. 
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3.  The  output  analysis/diagnosis  program  is  run  to  summarize 
and  present  the  results  of  the  analysis. 

OUTPUT  ANALYSIS /DIAGNOSIS  PROGRAM 

The  output  analysis/diagnosis  program  performs  the  following 

actions : 

1.  The  selected  intermediate  calculations  of  the  Areal 

Module  of  the  NRMM  for  a  single  terrain  unit  are 
transferred  to  the  diagnosis  subroutine. 

2.  For  each  slope  condition  (up,  level,  and  down)  the  output 

is  first  checked  to  determine  whether  the  terrain  unit 
is  a  GO  or  NO-GO  patch.  If  GO,  the  factor  which  limits 
the  speed  is  determined.  Otherwise  the  reasons  for 
the  NO-GO  is  determined.  In  either  case  a  code  is 
assigned  for  the  terrain  unit  and  slope  condition. 

3.  Steps  1  and  2  are  repeated  until  the  data  for  all  terrain 

units  have  been  read  and  analyzed.  An  output  file 
containing  the  terrain  unit  number,  speed-made-good, 
up,  level  and  down  speeds  and  diagnosis  codes,  and 
terrain  unit  area  for  each  terrain  unit  is  written. 

4.  The  area  (factor  area)  in  which  each  of  the  limiting 

factors  is  the  controlling  factor  is  determined  together 
with  the  average  speed  in  that  portion  of  the  total 
area  (for  the  GO  factors) .  More  precisely,  the  factor 
area  is  the  sum  of  the  areas  of  all  terrain  units  for 
which  the  factor  is  the  controlling  factor.  This  summary 
is  written  out.  (See  Appendix  A  -  Table  Al) . 
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5.  The  "Speed  Profile"  is  generated  for  each  slope  condition, 

i.e.,  the  terrain  units  are  sorted  into  decreasing  order 
of  speeds  and  the  cumulative  sum  of  the  area  and  average 
speed  are  computed.  The  speed  profile  data  are  written 
out  with  the  limit  code  for  each  terrain  unit  (generated 
in  Step  2) .  (See  Appendix  A  -  Table  A2) 

6.  The  speed  profiles  are  plotted.  On  the  plots  of  maximum 

speed  versus  percent  area,  different  symbols  are  used 
for  each  of  the  limiting  condition  codes.  (See  Appendix 
A  -  Figure  Al) 

Steps  1-3  are  most  efficiently  performed  during  a  run  of 
the  NRMM,  by  a  diagnosis  subroutine,  whereas  Steps  4-6  are  now 
performed  by  a  separate  program.  The  portion  of  this  program  which 
effects  Step  6  is  somewhat  non-portable  as  it  contains  calls  to  a 
plotter  routine,  of  necessity  device  and  system  dependent.  The 
remainder  is  standard  FORTRAN  and  the  plotting  section  is  quite 
simple.  Furthermore,  Step  6  only  presents  the  data  already 
available  from  Step  5  in  a  different  way.  However,  it  is  felt  that 
the  presentation  of  this  data  graphically  does  provide  a  good  way 
to  handle  the  data.  The  symbols  used  in  the  graphical  presentation 
on  the  plotter  are  an  arbitrary  set  available  conveniently  on  one 
system.  An  even  better  choice  would  be  a  color  coding  which  would 
be  easily  implemented  on  a  color-graphics  terminal,  but  this  presents 
difficulties  in  obtaining  hard-copy  output. 

The  output  of  the  analysis  program  consists  of: 

1.  A  list  of  any  terrain  units  for  which  the  program  was 
unable  to  determine  the  controlling  factor. 
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2.  The  summary  which  lists  the  factor  area  for  each  limiting 

factor. 

3.  The  list  of  all  terrain  units  in  decreasing  order  of 

speed-made-good  with  the  controlling  factor  code. 

4.  Speed  profile  plots  including  symbols  designating  the 

various  limiting  factors. 

In  the  case  of  a  terrain  unit  which  is  "GO"  for  the  vehicle 
and  slope  condition  (i.e.,  the  NRMM  predicts  a  non-zero  speed-made- 
good)  ,  the  analysis  is  based  on  various  sets  of  factors  which,  in 
combination,  give  rise  to  limiting  speeds,  one  for  each  set  of 
factors.  In  NRMM,  after  an  initial  screening,  candidate  speeds  are 
established  and  then  reduced  for  interactions  of  obstacles  and 
vegetation  with  the  sets  of  factors  first  addressed  (soil,  slope, 
ride,  tire  damage,  etc.).  Since  each  of  the  candidate  speeds  is 
an  upper  performance  limit,  the  final  speed  prediction  is  the  lowest 
of  the  various  limiting  speeds  (in  the  order  of  their  computation 
in  the  NRMM) .  In  the  output  Analysis/Diagnosis  Program,  the  diagnosis 
stops  when  one  of  these  candidate  speeds  is  equaled  (or  exceeded 
to  account  for  rounding)  and  the  corresponding  code  assigned.  The 
candidate  speeds  and  the  limiting  factors  to  which  they  correspond 
are  listed  in  Table  I. 

In  the  case  of  a  NO-GO  terrain  unit,  a  more  varied  collection 
of  variables  is  checked  to  determine  the  reason  for  the  NO-GO 
prediction.  Again,  these  are  assessed  in  the  order  of  their 
computation  and  evaluation  in  the  NRMM.  As  a  programming  convenience, 
the  values  assigned  to  the  limiting  factor  code  arc  negative  for 
the  NO-GO  analysis.  The  variables  used,  the  comparison  made,  and 
the  corresponding  reasons  are  listed  in  Table  II. 
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TABLE  I 

.  LIMITING  FACTOR  ANALYSIS  AND  CODE  ASSIGNMENTS 


SPEED 

CODE 

LIMITING  FACTOR 

VRIDE 

1 

Driver  tolerance  to  ride  over  rough 
terrain 

VTIRE 

2 

Tire  destruction  (applies  to  wheeled 
vehicles  only) 

. 99*VSOIL 

3 

Power  available  versus  resistance  due 
to  soil,  slope  and  overridden  vegetation 

VELV 

4 

Braking  available  relative  to  visibility 
restrictions 

V AVOID 

5 

Maneuvering  around  obstacles  and  overriding 
small  vegetation 

VBO 

6 

Maneuvering  around  vegetation  and  between 
obstacles 

VOLA 

7 

Obstacle  impact 

VOVER 

8 

Obstacle  and  small  vegetation  o\  ;rride 

WALK 

9 

Driver  prudence  in  vegetation  override- 

In  all  cases  where  the  variable  is  an  array  dependent  on  vegetation 
class,  the  value  used  for  the  comparison  is  that  of  the  vegetation 
class  used  to  compute  the  output  speed. 
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TABLE  II 

NO-GO  CAUSE  ANALYSIS  AND  CODE  ASSIGNMENTS 


VARIABLE  AND 

DECISION 

CODE 

NO-GO  CAUSE 

TBF  <  0 

-1 

Inability  to  brake 

VSOIL  £  0 

-2 

Soil  and  slope 

NEVERO  =  3 

-3 

Obstacle  interference 

NEVERO  =  1 

-4 

Belly  hangup  on  obstacles 

VBO  <  0 

-5 

Vegetation  too  dense  to  be 
avoided  and  too  large  to  allow 
override 

VXT  •-  d 

-6 

Tractive  force  needed  to  over- 
cide  obstacles  not  available 

The  value  vsoii,  used  here  is  that  without  vegetation.  The  values 
o:  VBO  and  VXT  are  those  of  the  vegetation  class  used  for  the  cor¬ 
respond  dig  computations  in  the  NRMM. 
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The  output  diagnosis  has  been  applied  to  runs  of  the  NRMM 
for  both  wheeled  and  tracked  vehicles  over  several  terrain  files, 
both  artificial  and  real.  All  of  the  terrain  units  have  been 
diagnosed.  (The  program  outputs  a  code  of  +99  to  designate  terrain 
units  for  which  the  limiting  factor  was  not  determined) .  The  analysis 
program  is  listed  in  Appendix  B.  The  additions  to  the  NRMM  including 
the  diagnosis  subroutine  required  are  presented  in  Appendix  C. 
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APPENDIX  A 


SAMPLE  OUTPUTS  OF  ANALYSIS/DIAGNOSIS  PROGRAM 


I 


UP  SLOPE  _  LEVEL  .  . PQUN_SLOPE _ 

LIMITING  MO. OF  FACTOR  ‘AVERAGE  NOTOF  FACTOR  AVERAGE  NO. OF  FACTOR  AVERAGE 
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TABLE  Al.  SUM'IARY  OUTPUT  FROM  OUTPUT  ANALYSIS/DIAGNOSIS  PROGRAM 


HP  Si  OF E  I  EVEL  .  _  . ItQU N  SLOPE _ _ 

"terrain' UNIT  CUMULATIVE  TERRAIN* UNIT  CUMLILATIUE  TERRAIn'uNIT  CUMULATIUE 

MAX  LIMIT  AREA  AVG  RCT  NO  MAX  LIMIT  AREA  AVG  F'CT  NO  MAX  LIMIT  AREA  AVG  F’CT 
SPEED  CODE  SPEED  AREA  SPEED  CODE  SPEED  AREA  SPEED  CODE  SPEED  AREA 
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TABLE  A  2 .  INDIVIDUAL  TERRAIN  L7IIT  OUTPUT  FROM  OUTPUT  ANALYSIS/DIAGNOSIS  PROGRAM 
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FIGURE  Al.  MAXIMUM  SPEED  UP  SLOPE  VS.  PERCENT 
OF  TOTAL  AREA  WITH  LIMIT  CODES 
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OF  TOTAL  AREA  WITH  LIMIT  CODES 
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FIGURE  A4.  MAXIMUM  SPEED  VS.  PERCENT  OF  TOTAL  AREA 

WITH  LIMIT  CODES  FOR  ALL  SLOPE  CONDITIONS 
(Original  in  Color) 
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APPENDIX  B 


LISTING  OF  OUTPUT  ANALYSIS  PROGRAM 


R-2183 


PROGRAM  PP’iRMM 


OUTPUT 


IALYSIS 


PROGRAM 


?:rmm  AREAL  MODULI 


DIMENSION  in(5£:S,3),V2(508,3),Al(530,3),AA(3), 

1  VlU(50  0),VlL(5DU),Vlu(5i?U),V2U(5U0),V2L(58n),V2D(5l33), 

2  PA,J(5ili*),PAL(5PI0),PAD(500),ICQDS(5?0,3),Nl(50'2,3) 
DIMENSION  NTGU(15,3),ATGO(15,3),VTGO(15,3),NT.NOGQ(13,3), 

1  ATNOGO  (  1 0, 3  ) ,  AGO (  3  ) , NGO  (3  )  ,  V  GO ( 3  ) ,  N'NOG 0( 3 )  ,  ANQGD  ( 3  ), 

2  LIflIT(15),NREAS(l0) 
nO'J3LE  PRECISION  L I M I T  ,  N  R  E  A  S 
EQUIVALENCE  (  7 1  (  1,  1 ),  V  l'J  (  1  )  ) 

EQUIVALENCE  ( V  1  ( 1 ,2  ), V 1L ( 1 ) ) 

EQUIVALENCE  ( V 1 ( 1, 3), VlD ( 1 ) ) 

EQUIVALENCE  ( A 1 ( 1 , 1 ) , P AU ( I ) ) 

EQUIVALENCE  ( A  1  ( 1 , 2 ) , P AL ( 1 ) ) 

EQUIVALENCE  ( Al ( 1 , 3 ) , P AJ ( 1 ) ) 

EQUIVALENCE  ( V2( 1, 1 ), V2U ( 1 ) ) 

EQUIVALENCE  (V 2(1,2), V2L(1)) 

E  QU I VALENCE  ( V 2 ( 1 , 3 ) , V  2D ( 1 ) ) 

DATA  (LIMI7(I),1=1,9)/4HRIDE,13IITIRE  C0NST,9HPQWER/RES, 

1 CHVI S I BILITY,  cliMAMSU  VER,  8!i'U NEU VER,  10HOBS  IMPACT, 9H03S  FORCE, 
1  JHDRVR.PRD.NC/ 


1000 

C 


DATA  (NREAS(I),I=l,o)/10HNQ  BRAKING, l/liSO IL/ SLOPE, 9H0B S  INTRF, 
1 0'siBELL  Y  HANG,  ID’!  VE  GET  AT  I  ON,  9  HOBS  FORCE/ 

ATOT=C . 

NPATC!i=0 
L 1 N=10 
LOUT=  22 

NP  ATCH=  N'P  ATCP  +  1 

INPUT  OF  SPEEDS  AND  CODES  FROM  DIAGNOSIS  SUBROUTINE 
READ  ( L I N  ,  £  N  D  =  l  v)  3  J  )  N  T  U  ,  1 7  U  T  ,  7.  M  ?  H 1 ,  (  V  1  (  N  ?  A  T  C  H  ,  J( )  ,  K  =  1 ,  3 ) , 

ZMPHS,  Z.«tPH6,2MPH7,ZM?Ha, GRADE,  AREA,  (ICODS  (KPATCH/K),  K=l,  3) 
FORMAT! I7,F10.4,3(F13.4,I4),Fie.4) 

DO  15  3  K= 1 , 3 

Al(. '.PATCH, K)  =  AREA 
N1  (  MPATCH/K  )  =  NTU 

CONTINUE 
ATOr=ATO?fAREA 
GUTO  133 

CONTINUE  - 

INITI ALIZA7I0II  \T- 

NP  ATCH  =  N?  ATCH- 1  *  * 

DO  1320  K  =  1 , 3 

DO  101/  L  =  1 ,15 


DO  101/  L=1 ,15 

NTGU(L,K)=0 
ATGO(L,K)=Z. 
VTGu ( L, K ) =U . 

CONTINUE 
DO  1015  L  =  1 , 10 

NT’.’OG  0 ( L, K ) -  ’ 
AT  1 0  (  L  ,  K  )  =  3  . 

CONTINUE 
NG0(K)=3 
N  N  0  G  0  (  K  )  =  J 
AGO(K 

A  N  0  G  0  (  K  )  =  0  . 

VG3(K )-A. 


u 

0 

\% 
O  yA 


1/20 


CONTINUE 


R-2183 


1050 

1052 


S  U  M  M  A  R  I  Z  A  ?  I  0  N 

Da  1  ] a u  i: = l , 3 

DO  13  53  1  =  1 ,  NP  ATCii 

L=ICOLS(I,X) 

IF(L.LT.  J)Gl)TJ  104? 

IF(  (L.  £2. 0)  .OR.  (L. GT. 15 )  )GGTO  1 25 *3 <5 
'/Tf;a(L/K)=:n,Ga(L/x)+i 
AT53(L/K)=AT51(L/K)+aI( I ,  K  ) 
VTOO(L,iO  =  VTCJ(L,K)+Vl(I,.\)*Al(I,X) 
!1GG(X)=NGU(X)+1 
AG0(X)=aG3(X)*A1(£,K) 

V30(K  )=VG0(X)*V1<I,K)*A1(I,K) 

GOTO  1060 

NO  GO 

IF(L.LT.-ITS)  GOTO  1053 
L=-L 

NTNi)Ga<L,X)=NTfJOGO(L,K)+l 
N’  N  0  G  0  (  K )  =  N  N  0  G  Q  (  K  )  + 1 
HT'!OGOCL/X)=AT*.T!GrKL/K)  +  Aia/:<) 
AN030(K)=AN0G3(K)+A1  (I,K) 

GOTO  1063 

ERROR  IN  ENCODING  OR  DIMENSIONS 

W  3 1  T  K  (  6  ,  1  3  5  2 )  L,N1(I,X),V1(1,K),K 
FORMAT (  '  ICUDF=', II,  '  NUT  RECOGNIZED'/ 

'  NTU, SPEED, SLOPE: ', 1 6, FI 3 . 2, 1 3/ ) 

CONTINUE 


DO  137?  L=l,15 

I E (  \TG3<L,K).SQ.0. )GOTO  1070 
V TGO(  L,  X  )  =  VTGO (L,  X  ) /  ATGO  (L , K  ) 

CONTINUE 

I F ( A G 0 ( K ) . S R . 3 .  ) GOT D  1100 
VGO( X ) =VGO( X ) / AGG( X ) 


1130  CONTINUE 

C  SUMMARY  OUTPUT 


WRITE (6/ 111?) 

1110  FORMAT(///33X,SHUP  SLOPE,  2U  X, 5LLE7EL, 1 9X, 1  DiiDOWN  SLOPS,  / 

1  1H  +  ,22X,3(4X,22  5I _ )/ 

2  2X/8HLI  !ITru;,13X,  3  (4  X,  5l:NO .  OF,  2X,  6H FACTOR ,  2 X, 

3  751  AVERAGE) / 3 X, 6 OF  ACTOR, 1 3X,  3(5  X,  5HT  U'S, 3X,4UAREA, 4X, 

4  5  H  S  P  E  E  D  )  / / / 1 7  0  GO  TERRAIN  UNITS/) 

DO  1140  L=1,V 

WRITE(6,113.;t)L,LT‘lIT(L),(LTGO(L,X),ATGQ(L,K),VTGO(L,K),K  =  l,3) 
113  0  FORMAT! IX, I  3, 1 X, A  1 0, 7 X, 3( II 3, 2F 8.2)) 

1143  CONTINUE 

:■!  k  I  T  E  (  6 , 1 1 5  F  )  ,  ( *  i  G  0  ( .<  )  ,  A  G  0  ( X  ) ,  7  G  0  {  X  )  ,  K  =  1 , 3 ) 

115  0  FORMAT (/IX, 5HTGTAL, 1 6 X, 3( 1 1 ?, 2FJ. 2  )  ) 


V  F<  I T  E  (  6 , 1 1 6  3  ) 

FORMAT!//'  LOGO  TERRAIN  UNITS'/) 


DO  llJG  L= 1 , 6 
LL=-L 

WR I TF ( 6, 1 1 7U )LL, NKF A 3 ( L ) , { N  TNOGG ( L , X  ), AT NOGO ( L, X  ), X  =  1 , 3  ) 
1170  FOR.XAT(iX,I3,lX,Alf),7X,3(llvi,Fe.2,aX)) 

1130  CONTINUE 

KR I TS( 6, 1 1 9  0 ) , ( NNOGO( X ) , A  NOGO( K ) , X  =  1 , 3 ) 

119  0  FORMAT (/IX, 5ETliTAL,16X,3( Il0,r3.2, 3X  )///  ) 

C 

B-2 
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□  F  SPEEDS 

N  L 1  =  N  P  A  T  C  H  - 1 
UJ  13  3  2  E  =  t  ,  3 

AT'-  ■'P=.‘ 
NTF.WP  =  3 
VT  EPP=3 
DO  13  32 


132^ 

1323 


IC  =  1, NL1 
ICP=IC+i 

no  1322  il=icp, npatcu 
1  F  (  » 1  (  IL/lO.La 
VTS/IP  =  V1(IL,K) 

A T E  1P  =  A 1  ( I L,  K ) 

UTEMP  =  N 1 ( IL, K ) 

I  T E  MP  =  I  C 3D  E  (  IL ,  X  ) 
V1(IL,X)=V1(IC,N) 
A1(IL,K)=A1(IC,!0 

M 1  (  IL,K)-N1(IC/X) 
ICr,DE(IL,<)=ICODE(IC,X> 
Vl(  IC,X)=VTEM? 

A  l  (  I C  ,  K  )  =  A  T  2  M  P 
M(  IC,.\)=N?2MP 
I C 1 D :•:  ( IC,X)=ITEMP 

CONTINUE 


71 ( IC,K ))GG  TO  1320 


% 

►— 

1/1 


CONTINUE 
DO  135  3  X  =  1 , 3 

VS  II ’■'=2. 
A  CU  •■'=2. 
33  1352 


t 

H 

3s 


8 


17  w 
£  *•» 
A  >  O 


1  =  1,  '.'PATCH 
A  C  U  M  =  A C  U  R  ♦  A 1  (  L,  K) 

V SU M=  V  SUM*  V 1 ( I , K ) *  A1 ( I / X ) 
41(I,K)=(ACUM/AT3T)*12a. 

I F (  VI  (  I,  X  ) .  LE .  3  . ) 30 TO  135Z 
V 2 (  1 , X  )  =V SUM/ AC HM 


1353 

1332 


1 

2 

3 

■1 


WRITE (22, 1330) 

FOPWAT  ( / 1 6X,  cril'P  SLOPE,  29X,  5HLE  VEL,  27X,  13iiD0W!l  SLOPE,  / 

1H  +  ,2X,3(2X,33H _ /3X)/ 

37,  3  (5X,  12::  TERR  A  IN  t!  NIT,  4X, 1 3HC UM UL A T I  V E, 3 X  )/ 

4X/3C3H  NO,  2X,3H  I A  X,  2X,  5liLI  !•!  IT,  1  X,  4H  ARE  A,  2  X,  JH  A  VS,  3  X, 

3**PCT,-*iX)/3X/3(5X,5MSPESD,2X4HCaD£,5X/bKS?SED/2X/4i!fcRcA/3X) ) 


r* 

DO  1410  I  =  l,NPATC;i 

D 3  1425  K=l,3 

\ 

DIFF=A1(I,K) 

£  F  (  l.  GT.  I  )  D  I  FF  =  D  I  FF- 
AA(X)  =  !)IFF*A?0T/133. 

1435 

CONTINUE 

W R I T E  (  2 2 ,  1  1 »’ C  )  ( N’  1  (I,K),V1(I, 
1  V2( I,K),A1( I,K),K=1 ,3) 

I  43  3 

FGRMAT(3(2X,I5,F6.1,1X,I3,F6.2,F6.1, 

1413 

CONTINUE 

- 

15  32 

T  7  P i  15  3  2 

FJ7MAT( *  ORAPLIC  OUTPUT  DESIRED?  7 

1512 

ACCEPT  1513,10 

FORNAX ( A 1 ) 

I  F  ( I  3  .  N  3 .  '  Y  *  )  5  T  0  P 

- 

C 

C 

A  M  A  X  =  A " A  X ! (VI (1,1), V 1(1, 2), V 1(1, 3)) 

SPEED  PROFILE  PLOTS 

c 

B-3 

.3  •  /  ■? 


J 


."JUS  ?lu2  LS  BJI5I 

mom  <xsn  fywiisKaC  to  bw 


2J13 

C 

c 

c 
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CALL  C  5  CL  (  A  M  A  X  ,  Y  v  A  X  ) 

CALL  GR  AF(  l  ,  .) .  ,  J  .  ,P  A’!,  VI  11,  NPATCi!,3  .  ,  l?  A  .  ,  V  . 

1  'PERCENT  TOTAL  AREA *,1 o, 'MAXIMUM  SPEED  I 

2  'MAXIMUS*.  SPEED  PROFILE,  J? 25,  1 ,  S,  3,  1  ,  : 
CALL  POINTS!  1,  NPATCii,71,  Al,  IC3DX) 

CALL  GRAF(-l) 

CALL  >  ?  A  F  (  D  .  ,  i P  ALT,  V  2M,  V  PATCH,  3  .  ,  l  D  3 .  ,  0  .  ,  ’ 
1  TL’TAL  ARM  A',  Id,  'MEAN  SPEED  (  "P'l )  1  6,  'A  if! 

221,1,5,  H, 1,0,  *PLT2',i)) 

CALL  J.r:  A  F  ( 1 ,  ,  3  .  ,  PAL  ,  Y1 1.,  N  ?  A  TO  H,  3  .  ,  1  33  .  ,  0  , 

1  'PERCENT  TOTAL  AREA ', Id,  'MAXIMUM  SPEED  I 

2  'MAXIMUM  SPEED  PRCFILE,  LF  V EL ',  28, 1 , 5,  .5, 
CALL  P0INTS(2, NPATCH, V1,A1,  I  CODE) 

CALL  GRAF(-l) 

CALL  GPAF (!.,?.,?  AL,  7  2L,  f.'PATCK,  fl.  ,  1 ,10.  ,  3  .  ,  • 
1  TOTAL  AREA',  Id,  'MEAN  SPEED  (:!?!!)',  1 5  ,' A  V  i 
221,1,5,0,1,0,  'PLT4',U) 

CALL  OR AF( 1 , , 0 . ,P  AD, VI D, NP ATCM, 2 • , 1 33 . , 0 , 

1  'PERCENT  TOTAL  AREA  ', 1 3,  'MAX  I UU“  SPEED  • 

2  '  "  A  X  I  M  U  M  SPEED  PROFILE,  DOWN  ',  2  7,  1,5,  i), 
CALL  POINTS  (3, NPATCH,  VI,  Al ,  I  CODE) 

CALL  GRA?(-t) 

CALL  CPkF  (  J  .  ,  ■).  ,?AC,V  2  D,  N?A  PC:!,  A.  ,  1  0  .  ,  ,3 . ,  3 
l  TOTAL  AREA', Id,  'MEAN  SPEED  (  •!!:’)  ' ,  IS,  'AVI 
221, 1,5, 8,1, (I,  'PLTM',0) 

FORMAT (  '  PEN  ',12,  'OSLO,  ESTER  PEN  ,;  ' ) 
FORMAT!  I) 

COMBINED  PLOTS  OF  ljP,LKVF.L  ■>  DOWN  PROFILES 


.,YMAX,10.,12., 
(MPH)  ',10, 

3,  VLT1  ',0) 


7 MAX, 13. ,13., ' PERCENT 
ERASE  SPEED  PROFILE', 

‘.,Y MAX, 13. ,13., 

( MP H )  ',10, 

!,!,«,  'PLT3', 3) 


'/MAX,  13.  ,13.,  'PERCENT 
ERASE  SPEED  PROFILE', 

!.,YMAX,1  3.  ,13.  , 
(MPa)', 19, 

1,2,  'PLTS',0) 


Y  -i A X  ,  1  3 .  ,1.1.  ,  'PERCENT 
.’ERASE  SPEED  PROFILE', 


CALL 

1  'i 

2  '1 
CALL 
CHAN’ 
I  P  =  N ! 
TYPE 

acce: 

I  P  2  = ! 

CALL 

CALL 

CHAN. 

I  P  =  N : 

TYPE 

ACC  E  ■ 

I  P  2  -  ’ 

CALL 

CALL 

CALL 

CALL 

1  ': 

2  '. 
CHAN, 
I  P  =  N 
TYPE 
ACCEi 
I  P  2  = 
CALL 
C  u  A  N 1 
IP  =  N 


GRAF  (1,0.,  3.,  PAN,  V1U,:.’  PATCH,  3.,  1  C  0  .  ,  ,  YM  A  X,  1  J.  ,  1 0.  , 

1 E R. C E ?i T  TOTAL  AREA', IB,  'MAXIMUM  SPEED  (MPH)  ',19, 
•AXIMUM  SPEED  PROFILE*, 21, 1,1,3. ,1,0, 'PLT7  ',1.  ) 

POINTS  (  1 , NPATCH,  VI,  A  1,1  CODE) 

'E  PEN 
•WPEN(3) 

2C  03,  I!’ 

’T  20 10, IP 

:evpem( ip) 

GRAF  (  0  ,  P  AL,  7  1L  ,  NP  ATCH,  1,7.) 

P  0 INTS ( 2 , NP ATCH, V 1 , A  1 , I CODE ) 


:wpen(  0  ) 

2300,  IP 
1  T  2310,1? 

IFWPFN  (IE) 

G?  AF ( 0 , PAD, V  ID, NP  A  TCN, 1,0) 

POINTS (3, NPATCH, VI, Al, l CODE) 

SR  AF (- l  ) 

GPAF(1,C.,  ,PAU,  V2U,  NPA  TCii,  3.  ,  1 Ufi .  ,  3 .  ,  YM  A X,  1  3.  ,  1  (3.  , 
'ERCENT  TOTAL  ARE  A  *,  1  d,  AN  SPEED  (MI-’M)',  15, 

OVERAGE  SPEED  PROF  ILE ',21,1,1,3.  ,1,2,  'PLT3  '  ,  1 .  ) 

;e  r  i 
:wf  o) 

2  3  J  -t ,  I P 

1  r  2/10,1? 

;ewpfn(  id 

GRAF (0, PAL, V2L,NP ATCH, 1,2) 

;e  pen 

:  WREN (3) 

B-.J 
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TVPE  220  1/  IP 
ACCEPT  22  IS, IP 

i?2=?iEv;psN(  ip) 

CALL  GRAF(0,PAD, V2J,NPA7CL,  1,3) 

CALL  GRAF  C  —  1 ) 

STOP 
?.  i  t  D 

S  J3R0UTINE  CSCL(A.M1X,  XMAX) 

scaling  subroutine  Fa1?  plot  axes 

I  F(  AUAX.LS.  '■!.  )  GO  TO  230 
FAC  =  1. 

AL5  =  ALGG1 3 ( 5 . ) 

AL  =  A  L  0  G 1 3  (  A  M  A  X ) 

EXP  =  AIST(AL) 

IF(2XP. EQ. AL)GO  TO  22 
IFUMAX.LT.  1.  )EXP=  EXP-1. 

R  =  AL-SXP 
FAC  =  13. 

IF(R.LE.AL5)FAC  =  5. 

XMAX  =  FAC  *  (13.  **EXP ) 

WRIT- (21,132) AMAX/AL, EXP, R, FAC, XMAX 
FORMAT ( IX, 63) 

RETURN' 

XMAX  =  3. 

RETURN 

END 

SU?HOJTiXS  POli’73  {K,\?ATCH,  VI  ,P  A,  l CODS) 
SUBROUTINE  FOR  OUTPUT  OF  COULD  DATA  POINTS  ON  SPE: 
DIMENSION  V 1  (  5  20, 3) , P  A( 52  3, 3 ) , I  CODE ( 500 , 3 ) 
DIMENSION  P  P  (  5  0  0  ) ,  7  P  (  5  0  3  ) 

00  2233  J=I,S 
NTU  =  3 

DO  2133  I  = 1 , NP ATCH 
I F ( I  CODE ( I , X ) . NE . J ) GO  TO  2133 
NT:J  =  NT-J  +1 
VP(NTU)  =  V 1 ( I , K ) 

PP(NTU)  =  PA(I,X) 

CONTINUE 

I F  (  N’TU . 30 . 0 ) GO  TO  2200 
CALL  GR AF( 0, PP, VP, NTU, -1, J ) 

CONTINUE 


PROFILE  PLOTS 


NO  GTS 

DO  2533  J=l, 6 
NTU  =  0 

DO  2433  I  = 1 , NPATCH 
IF(ICODE(I,X).NE.-J)GU  TO  2430 
NTU  =  NTU  +1 
VP  (NTU)  =  Vl(l,.<) 

PPCJTU)  =  ?  A  ( I  ,  X  ) 

CONTINUE 

I r ( NTU . 23. 3 ) GO  TO  2501 

CALL  GRAF ( 0, PP, VP, MTU , -1 , I 5-J ) 

CONTINUE 

RETURN 


Vi 


i 


t 


R-2183 


APPENDIX  C 


CHANGES  TO  NRMM  TO  IMPLEMENT  DIAGNOSIS 

1.  Add  array  ICODE  (with  DIMENSION  3)  to  COMMON  block  DERIVE.  DERIVE 

appears  in  the  main  program  and  subroutines  AREAL,  ROAD, 
and  BUFF0. 

2.  Add  the  control  variable  KDIAG  to  the  COMMON  block  SCEN.  SCEN 

is  in  the  main  program  and  subroutines  SCN ,  VPP,  TERTL, 

AREAL,  ROAD  and  BUFF0. 

3.  Add  KDIAG  to  the  NAMELIST  CONTRL,  lines  SCN16-SCN25  of  subroutine 

SCN.  This  allows  KDIAG  to  be  read. 

4.  In  subroutine  SCN,  change  line  SCN-55  to: 

IF  (DETAIL. EQ. 3}  GO  TO  313 
Add  between  lines  SCN  65  and  SCN  66  the  following: 

313  CONTINUE 

KDIAG  =  1 
GO  TO  330 

These  changes  permit  diagnosis  to  be  specified  as  output  detail 
level  3. 

5.  In  subroutine  AREAL,  add  between  lines  AREAL-219  and  AREAL-220 

the  following: 

IF  (KDIAG. EQ.l)  CALL  DIAG  (AREA,  IDCODE,  IMAX,  ITUT,  NEVERO, 
NTU,  TBF,  VAVOID ,  VBO,  VELV,  VMAX,  VMAXl ,  VOLA,  VOVER,  VRID, 
VSOIL ,  VTIRE,  VXT) 

6.  Change  line  BUFF0  34  of  subroutine  BUFF0  i_o: 

ZMP1I7 ,  ZMPH8,  GRADE,  AREA,  ICODE. 
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This  addition  to  the  binary  output  file  written  to  LUN10  allows 
this  output  file  to  be  used  by  the  analysis  program. 

7.  The  diagnosis  subroutine  DIAG  must  be  included.  The  listing 
of  this  subroutine  follows. 
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S  U  8  R  0  U  T IV  E  UIAG(ARF»,lCnDS,iy.AX,t?VT, 5  V  "  3,  X  T  V,  TO  F, 

1  VAVOIQ/VBO/VELV, VMAX, VMAXl/VCLA,  VjVER,VRID,VSUIL 

2  VT1RE,VXT,V.»:ALK) 

DIAGNOSTIC  SUB ROUT  I  f.'E  rOk  NRMM  AREAL  MGOULE 

DIMENSION  ICOJE(3),IHAX(3),TBF(3),VaV(IIu(3,9),VjO(  3,9)/ 

1  VEL V ( 3  ), VM  AXl  (  3 ),  VO  V  'R  (  3,  9  ),  VGQ I L  (  3, 9  ) ,  V  XT  (  3, 9  ) 

IF(1T;JT.GE.11)GUT0  1333 

assignment  of  limiting  kactur 
do  x  - 1 , 3 

VC  =  VVAXl(K) 

IF(  VC.EO.  3.  )GO  TC  4  ).1 

HATCH  IS  SO 

IF(7C.LT.  VRIOGO  TC  213  \ 

.TIDE  LIMIT 

lCODE(K)=i 
GO  TO  603 

I  F  (  V  C  .  N  E  .  V  T  i  R  E  )  GO  TU  223 
TIRE  CONSTRUCTION  LIMIT 
ICODF (X)=2 
GO  TO  6-33 

IFI  7C.LT.  (•9r>*V30IL(;\/  IMAX(K))  )  )GG  TO  233 
POWER  VS.  SOIL, SLO°E, VEGETATION  RESISTANCE 
I  CODE ( K )  =  3 
GO  Ttl  603 

IFCVC.LT.  VEL7(K))CO  TO  243 
V1SI3LIT7  LIMIT 

I  CODE ( K )  =  4 
SO  TO  6J3 

IFCVC.LT.  VAV3IC(K,I  !AX(K)))SO  TO  253 
MANEUVER  ARO'JN'D  OBSTACLES  A  NO  VEGETATION 
I  CODE ( K ) =5 
GO  TO  60S 

IFCVC.LT. V3C(K,1MAX(K)))GJ  TO  263 
MANEUVER  AP.OUNL  VEGETATION 
ICODF (K)  =  6 
GO  TO  633 

IF( VC. LT. VCLA ) GO  TO  273 
G3STACLE  IV.^CT  LIMIT 
I  CODE (K )  =  7 
GO  TU  633 

IFCVC.LT.  VOVERIX,  IV  AX  (X  )  ) )  GOTO  283 
POWER  TQ  OVERRIDE  OBSTACLE 
ICODF (K)=3 
GOTO  603 

DRIVER  PRUDENCE  OVERRIDING  VEGETATION 
I-CVC.NE.V WALK ) GOTO  2  93 
I  CODE ( K  )  =  9 
GOTO  663 

LIMIT  NUT  DIAGNOSED 
lCODE(< )  =  9D 
GO  TO  633 
CONTI NUE 

PATCH  IS  NO  GO 


IF(ThFU)  .SL.U.k5)C,0  Til  413 
'JO  BRAKING  NO  GO 

ICODF ( K ) =- I 
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GO  TO 

41?  1F(  VS01L(:C,1).CT.3.0)C0  TO  423 

c  soil  i.  scoff  if-iosiLiLATio:; 

I  CODS  (K)=~2 
GO  TO  60 

42?  IFCiE VFRO.NL.  3)30  TO  43? 

C  C3S7ACLS  INTERFERENCE 

rCi!0f:(K)--3 
GO  TO  622 

432  IF (NF.VERQ.NE.  1  )GCTG  4 *12 

C  BELLY  HaNSU?  ON  OBSTACLES 

I  CODE (  K  )  =  -  1 
GOTO  623 

44?  l?(VhG(K,  IMAX{K)).GT.O.  )GGTO  40.5 

c  vegetation  to  go 

I  CODE ( X  )=-5 
GOTO  62.5 

45?  IF(  (  NSVERG.GT.  £? ) .  OR.  (  VXTCK,  I  MAX  (<  )  ).  GT.  3 . 5  JGOTO  463 

C  OBSTACLE  OVFRRIJE  FORCE  MOGO 

I CODE (K )=-6 
GOTO  633 

C  tIGGO  REASON  LOT  DIAGNOSED 

46  3  IC0DE(? )=-99 

603  CONTINUE 

C  WRITE!  21 , 61 2 ) NTU,  V!  AX, ( VMAX1 ( * ) , I  CODEC  N ) , 4=1 , 3 ) , ARE A 

61?  FORMAT  (IX,  I6,F14.  4,3(r  12.4,  I4),F13.4) 

RETURN 

1233  WRITE(6,li3n) 

1210  FORMAT  ( 42;!  POAD  TSPRAIN  UNIT  DIAGNOSIS  NOT  AVAILABLE) 

END 


C-4 


